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ABSTRACT 



> 

o 
o 

■ Here we report on the first XMM-Newton observation of the Anomalous X-ray 

O ! Pulsar IRXS J170849.0-400910. The source was observed in 2003 August and was 

. found at a flux level a factor of about two lower than previous observations. Moreover, a 

' significant spectral evolution appears to be present, the source exhibiting a much softer 

fH \ spectrum than in the past. Comparison of the present properties of IRXS J170849.0- 

• 400910 with those from archival data shows a clear correlation between the X-ray flux 

and the spectral hardness. In particular, the flux and the spectral hardness reached a 
maximum level close to the two glitches the source experienced in 1999 and in 2001, 
and successively decreased. Although the excellent XMM-Newton spectral resolution 
should in principle allow us to detect the absorption line reported in a phase-resolved 
spectrum with BeppoSAX. , and interpreted as a cyclotron feature, we found no ab- 
sorption features, neither in the phase averaged spectrum nor in the phase resolved 
spectra. We discuss in detail both the possibilities that the feature in the BeppoSAX 
data may have resulted from a spurious detection or that it is real and intrinsically 
' variable. We then discuss a possible explanation for the glitches and for the soften- 

ing of the source emission which followed the flux decrease, in the framework of the 
magnetar model. 
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fields - stars: neutron - X-rays: stars 



1 INTRODUCTION 

Anomalous X-ray Pulsars (AXPs) and Soft 7-ray Repeaters 
(SGRs) are two peculiar groups of neutron stars (NSs) which 
stand apart from other known classes of X-ray sources. They 
are all radio-quiet, exhibit X-ray pulsations with spin pe- 
riods in the ~5-12 s range, a large spin-down rate {P « 
10"^'' — 10~^^ss~^) and a rather high X-ray luminosity 
{Lx ~ lO"^* — 10"^®ergs~^; for a recent review see Woods 
& Thompson 2004). The nature of the X-ray emission from 
these sources has been intriguing all along. In fact, for both 
AXPs and SGRs, the X-ray luminosity is too high to be pro- 
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duced by the loss of rotational energy of the neutron star 
alone, as for more common isolated radio pulsars. 

The magnetic fields of AXPs and SGRs, as estimated 
from the classical dipole braking formula B ~ 3.2 x 

10"\/P? G, are all above the critical magnetic field at 
which quantum effects become important, Bqed ~ 4.4 x 
10^^ G. At the same time, the lack of observational signa- 
tures of a companion strongly argues against an accretion 
powered binary system, favoring instead scenarios involv- 
ing isolated NSs. These findings led to the idea that these 
two classes of sources should be somehow linked together, 
and their X-ray emission related to their very high magnetic 
fields. 

At present, the model which is most successful in ex- 
plaining the peculiar observational properties of AXPs and 
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Figure 1. Compajrison of the folded lightcurves of IRXS J170849.0-400910 as a function of energy for the BeppoSAX (Rea et al. 2003, 
left) and XMM-Newton EPIC-PN (this paper, right) observations. Note that due to the time elapsed between the two observations, the 
pulse profiles are not phase connected. The solid lines indicate the minimum of the 0.5-1 keV in the two cases and is superimposed to 
the pictures in order to better show the shift in phase of the two sets of folded lightcurves at different energies. Units on the y-axis are 
counts per seconds with an arbitrary normalization. 



SGRs is the "magnetar" model. In this scenario AXPs and 
SGRs are thought to be isolated NSs endowed with ultra- 



high magnetic fields (B 



10' 



G) and their X-ray 



emission is powered by the magnetic field decay (Duncan & 
Thompson 1992; Thompson & Duncan 1993, 1995, 1996). 
Alternative scenarios, invoking accretion from a fossil disk 
remnant of the supernova explosion (van Paradijs, Taam 
& van den Heuvel 1995; Chatterjee, Hernquist & Narayan 
2000; Alpar 2001), encounter increasing difficulties in ex- 
plaining the data. 

While SGRs were first identified through intense, re- 
peated high-energy bursts observed in the X/7-ray band, 
AXPs were revealed thanks to their persistent omission in 
the soft X-rays. For more than two decades, AXPs were 
thought to be steady, soft X-ray emitters and their X-ray 
spectrum was well modeled as the superposition of a black- 
body at kT ~ 0.4 keV and a steep power-law with photon 
index P ~2-4. No counterparts had been detected at other 
wavelengths. Only in the past few years, mostly thanks to 
new generation observatories and to dedicated observational 
campaigns, our picture of AXPs has largely changed: 1) a 
few of them show definite long term X-ray flux variability 
on timescales of months (Woods et al. 2004; Mcrcghctti et 
al. 2004; Gavriil & Kaspi 2004), 2) large spectral variations 
with rotational phase were revealed (Israel et al. 2001; Rea 
et al. 2003; Tiengo et al. 2005), 3) bursting activity was de- 



tected (Kaspi et al. 2003; Gavriil, Kaspi & Woods 2002), 
4) glitches were discovered (Kaspi, Lackey & Chakrabarty 
2000, Dall'Osso et al. 2003; Kaspi & Gavriil 2003, Morii, 
Kawai & Shibazaki 2005), 5) high-energy tails requiring 
an additional spectral component were discovered (Kuiper, 
Hermsen & Mendez 2004) and 6) IR/Optical steady and 
variable emission was detected (Hullcman, van Kerkwijk & 
Kulkarni 2000; Kern & Martin 2002; Wang & Chakrabarty 
2002; Israel et al. 2002). 

IRXS J170849.0-400910 was discovered with ROSAT 
(Voges et al. 1996) and later on a ~ lis modulation was 
found in its X-ray flux with ASCA (Sugizaki et al. 1997). 
Early measurements suggested that it was a fairly stable ro- 
tator (Israel et al. 1999). However, in the last four years the 
source experienced two glitches, with different post-glitch 
recoveries (Kaspi, Lackey & Chakrabarty 2000, Dall'Osso 
et al. 2003, Kaspi & Gavriil 2003). Searches for optical/IR 
counterparts ruled out the presence of a massive companion 
(Israel et al. 1999), even though very recently an IR coun- 
terpart has been proposed (Israel et al. 2003; Safi-Harb & 
West 2005). A diffuse (~ 8') radio emission at 1.4 GHz was 
recently reported, possibly associated with the supernova 
remnant G346.5 0.1 (Gaensler et al. 2000). 

Pulse phase spectroscopy analysis of two 
BeppoSAX observations (Israel et al. 2001; Rea et al. 
2003) revealed i) a large spectral variability with the spin- 
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Figure 2. The 0.5-10 keV EPIC-PN phase-average spectrum of IRXS J170849.0-400910 . The residuals are referred to a fit with an 
absorbed blackbody plus a power-law. Systematic deviations are present in the residuals, most probably due to instrumental effects. 



phase, ii) a strong energy dependence of the pulse profile 
shape, and iii) shifts in the pulse phase between the low and 
the high energy profiles, ffigh variability of the pulse shape 
with energy is now detected at even higher energies, up to 
~220 keV, both with the HEXTE instrument on board of 
RXTF, and the INTEGRAL satellite (Kuiper et al. 2005, 
in preparation). 

By analyzing a 5eppoSAX observation taken in 2001 
(the longest pointing ever performed on this source), Rea 
et al. (2003) reported the presence of an absorption line at 
~ 8 keV in a phase-resolved spectrum. The line significance 
is ~ 4cr. Interpreting the feature as a cyclotron line due to 
resonant scattering yields a neutron star magnetic field of 
either 9.2 x 10^^ G or 1.6 x lO^'^ G, in the case of electron or 
proton scattering, respectively. 

In this paper we report on the first XMM- 
Newton observation of the AXP IRXS J170849.0-400910 . 
Timing and spectral data analysis are presented in §|5| The 
Discussion is divided in two parts: in Section 3 we discuss the 
results, readdress the problem of the line significance in the 
BeppoSAX observation, and briefiy report on the cyclotron 
line variabilities detected so far in other X-ray pulsars; Sec- 
tion 4 is instead a more theoretical part where we propose 
an interpretation of the results in the framework of the mag- 
netar scenario. Conclusions follow in §|S| 



2 OBSERVATION AND DATA ANALYSIS 

IRXS J170849.0-400910 was observed with XMM- 
Newton between 2003 August 28th and 29th, for ~ 50 ks. 
The XMM-Newton Observatory (Jansen et al. 2001) 
includes three 1500 cm'^ X-ray telescopes with an EPIC 
instrument on each focus, a Grating Spectrometer (RGS; 
der Herder et al. 2001) and an Optical Monitor (Mason 
et al. 2001). Two of the EPIC imaging spectrometers use 
MOS CCDs (Turner et al. 2001) and one uses a PN CCD 



(Striider et al. 2001). The MOS cameras were operated in 
Prime Partial Window Mode, while the PN camera was in 
Prime Small Window Mode, all with the medium optical 
photons blocking filter. The Optical Monitor was operating 
during the observation but the optical counterpart of 
IRXS J170849.0-400910 was too faint to be detected. Data 
were processed using SAS version 5.4.1. We employed 
the most updated calibration files available at the time 
the reduction was performed (June 2004). Standard data 
screening criteria were applied in the extraction of scientific 
products. Since a higher background affected the last ~10ks 
of the observation, we used only the data during intervals 
in which the count rate above 10 keV was less than 0.35 
counts s~^. The source events and spectra were extracted 
within a circular region of 27" centered on the peak of the 
point spread function of the source. This non standard 
radius was used because the source was located near the 
edge of the chip. The background was instead obtained from 
a source- free region of 27". When generating EPIC spectra 
we have included events with PATTERN 4 (i.e. single and 
double events). We have checked that spectra generated 
with only single events (i.e. PATTERN= 0) agreed (apart 
from normalization factors) with the spectra generated 
from single and double events. All of the EPIC spectra 
were rebinned, before fitting, using at least 20 counts per 
bin and not oversampling the resolution by more than a 
factor 3. Since both MOS and PN give consistent results 
as far as the timing analysis, and MOS data are slightly 
affected by pile-up, in the following we only report on PN 
data. We extracted first and second order RGSl and RGS2 
spectra for both source and background using the standard 
procedure reported in the XMM-Newton analysis manual 
but found no evidence for spectral features. 

Thanks to the high timing and spectral PN resolution 
we have been able to perform timing and spectral analysis, 
as well as Pulse Phase Spectroscopy. Results are reported in 
the following subsections. 
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Figure 3. Spectral parameters as a function of phase. The open circles represent the fits with A'^^f as a free parameter, while the dashed 
lines indicate the phase-resolved spectral parameters with Nfj fixed at the phase-average value. 



2.1 Timing Analysis 

In order to determine the spin period of IRXS J170849.0- 
400910 we first barycentered the events (using the SAS tool 
barycen) and performed a power spectrum (powspec) where 
we found a periodicity at ~lls (as expected from previ- 
ous period determinations) composed of a fundamental and 
a second harmonic. In order to obtain a precise estimate 
of the pulse period, we carried out a period search around 
this value (ef search). We obtained a best spin period of 
Ps = 11.00170 ± 0.00004 s. The uncertainty was determined 
by dividing the observation in 8 intervals and performing a 
linear fit to the phase determined in each interval (phase- 
fitting technique). 

The PN phase-folded light curve in different energy 
bands is shown in Fig0 (right panel), together with that ob- 
tained in the past by BeppoSAX . We found that the pulsed 
fraction of the X-ray signal (defined as the amplitude of the 
best-fitting sine wave divided by the, background corrected, 
constant level of the emission) is energy- dependent, and it 
varies from 39.0 ±0.5% in the 0.5-2.0 keV range to 29 ±1.5% 
in the 6.0-10.0 keV range. These values are consistent with 
those reported for the pre-glitches BeppoSAX observation 
(Israel et al. 2001) while both are larger than those re- 
ported for the post-glitches BeppoSAX observation (Rea et 
al. 2003). 

2.2 Spectral Analysis and Pulse Phase 
Spectroscopy 

We generated the PN phase-averaged spectrum and 10 
phase-resolved spectra. Phase zero was arbitrarily defined 



as the start of the observation (MJD = 52879.9115833) and 
the phase bins are equally spaced. The choice of the number 
of intervals was based on the comparison between the better 
quality of the A'MM-A''eu)ton phase-averaged spectrum with 
respect to the BeppoSAX spectrum (for which 6 phase in- 
tervals were used) and it was made prior to the analysis, 
in order to keep the number of trials to a minimum in case 
evidence for a cyclotron line was present. 

The phase-averaged spectrum in the range 0.5-10 keV 
is satisfactorily fitted by a model consisting of an absorbed 
blackbody plus a power law (xlo/ — 1-28 for 242 degree of 
freedom, with a 2% systematic error included; see Fig|2|l. 
The best fitting parameters are: hydrogen column density 
Nh = (1.48 ± 0.04) X 10^^ cm"^, blackbody temperature 
fcTbb = 0.456 ± 0.009 keV (with a blackbody radius of 
Rhh = 2.7 ± 0.2 km at a distance of 5kpc) and photon in- 
dex r = 2.83±o:gi (all errors are at the 90 % confidence 
level)^. The unabsorbed flux in the 0.5-10 keV range is 
9.1 X 10~^^ erg cm""^ s^^ corresponding to a luminosity of 
2.7 X 10^^ ergs~^ (again for a 5 kpc distance). In the 0.5-10 
keV band, the blackbody component accounts for ~ 14.9 % 
of the total unabsorbed flux. We found some deviations from 
the best fit model in the lower energy range (see Fig|2J , how- 
ever we believe these deviations are the result of mixing of 
spectral shapes, since those deviations are not present in 
the phase-resolved spectra and to calibration issues. Consis- 
tently with what has been done in the past for other AXPs, 



^ Note that using the XSPEC phabs model instead of the wabs 
used here, Nh = (1.36 ± 0.03) X 10^^ cm-^ while the other pa- 
rameters do not change significantly. 
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Figure 4. Residuals (tr) of the fits to tiie phase-resolved spectrum of IRXS Jf 70849. 0-4009f0 . Phase intervals are indicated in the 
lower-right of each panel. This figure shows the absence of any significant feature during the whole pulse period. 



we tried to fit the spectrum of IRXS J170849.0-400910 with 
two absorbed blackbodies but this resulted in a worse fit 
ixlof ~ 2.2). 

The absorbed blackbody plus power law model provides 
excellent fits for all the ten pulse phase- resolved spectra, ei- 
ther when leaving Nh free or when fixing it to the value 
obtained from the phase-average analysis. The resulting pa- 
rameters (both for Nh fixed and free) are shown in Fig. |3] 
We note that the lowest is observed when the blackbody 
fiux is lowest, which might suggest that a simple blackbody 
is not a very realistic description of the thermal component. 
As we tested, by replacing the single temperature blackbody 
by a disk blackbody gives a lower Xdof (reduced by up to 
0.06) in 9 out of the 10 phase-resolved spectra. Of course 
this does not imply much about a physical interpretation, 
but is suggestive of a thermal emission more complex than 
a single blackbody. 

Note that all the phase-average spectral parameters, 
but the power-law photon index F, are consistent (within 
1.5a errors) with the last 5eppoSAX observation. We then 
found here a clear softening in the spectrum correlated with 
a decrease of the soft X-ray flux by a factor of almost two. 

As we can see from Fig^] there is no significant indica- 
tion for an absorption line near 8.1 keV, as reported by Rea 
et al. (2003). In order to determine an upper limit to the 
strength of such a line, we concentrated on the phase inter- 
val which coincides best with that in which the absorption 
line was found in the 5eppoSAX data. From the re-analysis 
of the 5eppoSAX data, we noticed that the phase inter- 
val in which the absorption line was strongest is incorrectly 
stated in the published version of Rea et al. (2003) (see also 
S 13.111 . The correct value is 0.82-1.0 in their phase conven- 
tion (phase zero corresponding to minimum of the 0.1-2 keV 
pulse profile) which corresponds best to phases 0.5-0.6 and 
0.6-0.7 in our phase convention. Thus, it does not lie in the 
rising part of the pulse profile, but close to the minimum 
in the 0.1-2 keV range (which actually is the maximum in 



the 6-10 keV band; see Fig0. We then fitted the spectrum 
taken at these phases by adding to the best fitting contin- 
uum (with Nh fixed) a cyclabs model (Mihara et al. 1990; 
Makishima et al. 1990) as that used by Rea et al. (2003). We 
fixed both line energy and width to 8.1 keV and 0.2 keV, re- 
spectively, finding no significant improvement in the ■ We 
stepped through values of the dimensionless line depth while 
the continuum parameters were allowed to change, until a 
Ax^ of 2.706 (90% uncertainty region for one parameter of 
interest) was reached. This occurred at a depth of 0.13 and 
0.16 for the two phase intervals, and we consider the value 
0.15 as a 95% upper limit to the line depth. 

Hints for spectral lines could be found around ~ 7keV 
in some phase-resolved spectra and around 1.8 keV and 
2.8 keV in the phase-average spectrum (see Fig|l]and Fig|5J 
but in all cases the significance considering one single trial 
was less than 2a, then highly not statistically significant. We 
believe them most probably due to systematic instrumental 
deviations. 



3 DISCUSSION ON THE RESULTS 

In this paper we analyzed the first XMM- 
Newton observation of the AXP IRXS J170849.0-400910 . 
We found the source at a flux level two times lower than 
that measured in a previous BeppoSNK observation (Rea 
et al. 2003) and with a significantly softer spectrum (see 
FiglSl^ )• The strong phase dependence of the spectrum, the 
large pulse shape variability and the shifts in phase of the 



^ Note that the BeppoSNX. flux reported in this plot is slightly 
lower than one reported in Rea et al. 2003, this is due to a dif- 
ferent, more updated, arf matrix used here for the re-analysis of 
this data. 



6 Rea et al. 




3 4 5 6 7 8 

Energy (keV) 



Figure 5. Monte Carlo simulations. Cyclotron line depth as a 
function of energy for 10000 simulated spectra. 

pulse profile as a function of energy, previously detected, 
are confirmed by tlie present data. Comparing the spectral 
parameters obtained from this XMM-Newton ohserY&tion 
with those derived from the BeppoSAX observation, we 
find that Nh and blackbody temperature are consistent 
with being constant, while the blackbody radius has de- 
creased by about a factor ~2.4. At the same time, the flux 
decreased substantially (by a factor 2) and the power-law 
index increased from ~ 2.4 to ~ 2.8. 

The long-term evolution of the source flux and the spec- 
tral hardness are shown in Fig|S| where different observa- 
tions, spanning nearly flve years, have been collected The 
correlation of the photon index with the source luminosity is 
evident by comparing the two panels. The spectrum became 
progressively harder as the flux rose in correspondence of the 
two glitches and then softened as the luminosity dropped, 
following the glitch recovery. This is suggestive of a scenario 
in which the mechanism responsible for the glitches is also 
at the basis of the enhanced emission and of the spectral 
hardening. 

The similarity of the second glitch of IRXS J170849.0- 
400910 with the one discovered during the bursting activity 
of IE 2259+586 (Kaspi et al. 2003), after which a similar 
exponential recovery was seen, suggests that bursts likely 
occurred in IRXS J170849.0-400910 as weU, but the sparse 
observations did miss them, as already suggested by Kaspi 
& Gavriil (2003). Moreover, the spectral parameters and 
the flux changes after the recovery of the glitch strengthens 
this idea in comparison with what has been reported for the 
post-bursts fading of IE 2259+586 (Woods et al. 2004). 

3.1 The absorption line: Spurious vs. Variability 

As discussed above, we have found no evidence in this 
XMM-Newton observation for an absorption line at ~ 
8 keV neither in the phase-averaged nor in the phase re- 
solved spectra; the upper limit for the line depth is 0.15 
at 95% confidence level. If we compare this with the value 



found by Rea et al. (2003; 0.8±0.4 at 90% confidence level), 
there is only a very small chance that the two measurements 
are consistent. This leaves us with two options to be consid- 
ered: i) the BeppoSAX detection is a spurious result, or ii) 
the absorption line is variable. 

Before discussing the option of a variable absorption 
line (see JIJ, we undertook a careful re-analysis of the 
5eppoSAX data. This re-analysis resulted in the finding 
that the phases at which the absorption line was strongest 
were given incorrectly in the published version of Rea et 
al. (2003), as we noticed earlier: in particular, the line is 
strongest close to the pulse minimum in the 0.1-2 keV band 
(or the pulse maximum in the 6-10 keV band). We never- 
theless found that the reported estimate of the significance 
is sound and not much influenced by different choices in 
the background subtraction (annular regions, circular re- 
gions far from the source or using blank field files) or by 
different extraction regions for the source. The re-analysis 
of the BeppoSAX data made varying the extraction radius, 
the criterion for the background subtraction and the spectral 
binning factor, results in basically the same line properties, 
which strengthens our confidence in the robustness of the 
result. 

Using an F-test method and taking into account the six 
trials we made in the phase resolved spectra, we derive a con- 
fidence level for the absorption line of ~ 4(t. Note that even 
if we take into account all the possible energies at which the 
feature could lie in the LEGS plus MEGS energy range, the 
confidence level is still ^ 3.5o". We note that recent works 
pointed out that the F-test may be inappropriate in this cir- 
cumstances (see Protassov et al. 2002) , leading sometimes to 
incorrect estimates of the significance of a feature, althought 
it has been (and still is) widely used to test the significance 
of spectral lines. 

In order to further investigate on this issue, we ran a 
Monte Garlo simulation of 10000 spectra with only the con- 
tinuum model (parameters reported in Rea et al. 2003) and 
the same number of photons of the phase resolved spectrum 
which showed the line in the 2001 BeppoSAX observation. 
The results of the simulation is shown in In Fig|3 Each circle 
represents the line depth found in each of the 10000 spectra, 
here plotted as a function of the line energy. From this simu- 
lation we found 32 spectra with depth >0.8 in 10000 points. 
We can then reliably say that the probability of the line 
being a fiuctuation is <0.32%. In summary, we confirm the 
detection of the 8.1 keV absorption line in the BeppoSAX 
data made by Rea et al. (2003) at 99.68% confidence level. 

The interpretation of the absorption feature as a cy- 
clotron scattering line proposed by Rea et al. (2003) was 
based on the following criteria: 1) a Gaussian line gives a 
bad fit and does not reproduce the asymmetrical shape of 
the observed feature; 2) the best fitting model is the XSPEC 
cyclabs model; 3) the line strength is highly phase depen- 
dent; 4) no atomic edges or absorption lines are known to 
lie around 8.1 keV (at least without assuming ad hoc shifts 
possibly due to the high gravitational redshift or to Zeeman 
effects in such strong magnetic field); 5) the relation between 
the line energy and width agrees with that of cyclotron scat- 
tering features discovered in other classes of sources (see 
Fig|H]in Rea et al. 2003); 6) the magnetic field inferred from 
the line energy, either in the case of an electron or proton 
cyclotron resonance, is reasonably consistent with what is 
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expected for a normal neutron star (~ 10^'^ G) or for a mag- 
netar (~ 10^^ G), both being still open possibilities. Then, 
if this feature is real, all the above points hint toward the 
cyclotron nature of the absorption feature at 8.1 kev. 

Keeping always in mind the possibility that the absorp- 
tion line in the BeppoSAX spectrum might be due to sta- 
tistical fluctuations, in the following section we discuss a 
physical mechanisms which could be responsible for the ap- 
pearance of a transient cyclotron line in this source in the 
context of the magnetar scenario. 



and has been detected only in the spectrum of a number of 
quite energetic bursts. Moreover, recently a hint for a 4 keV 
feature has been found during bursts detected by XMM- 
Newton (Mcrcglietti et al. 2005). Despite intensive searches 
no evidence of any feature has ever been found in the per- 
sistent X-ray emission of this source. The feature appears to 
be almost certainly correlated with the onset of the bursting 
activity and with periods of increased X-ray luminosity and 
spectral hardening. 



3.2 Other cases of variability of cyclotron 
absorption features 

We conclude this section by summarizing the current star- 
tus of detections of variable resonant scattering cyclotron 
absorption features in X-ray pulsars, both in binary and in 
(probably) isolated systems. 

Gyclotron resonant scattering features were first de- 
tected in accreting X-ray pulsators. Up to now more than 
ten pulsars in binary systems axe known to exhibit an ab- 
sorption line in their spectra. In all these cases, there is 
unanimous consensus about an interpretation of the feature 
in terms of resonant scattering by electrons (see e.g. Heindl 
et al. 2004 for a recent review). Variability in the line energy 
has been observed in some of these sources, and changes ap- 
pear to correlate with the source flux. This behavior finds a 
quite natural explanation within the column accretion sce- 
nario, in which the height of the radiation shock depends on 
the luminosity. At larger fluxes the shock forms higher up 
in the column where the (dipole) fleld is weaker (Mihara, 
Makishima & Nagase 1997; Mihara, Makishima & Nagase 
2004). 

Cyclotron features have been reported so far in the spec- 
tra of some isolated neutron stars too: the soft 7-repeater 
SGR 1806-20, flve X-ray dim isolated neutron stars (XDINs) 
and the radio-quiet neutron star IE 1207.4-5209. With the 
exception of IE 1207.4-5209 (Bignami et al. 2003), the na- 
ture of the line (either electron or proton) is controversial 
as yet, and, at least for XDINs, a possible atomic origin can 
not be excluded. 

The similarity of XDINs pulsational periods (^ 3-10 s), 
together with the alleged high values of their magnetic field 
(- 10^^-10" G), makes comparison with IRXS J170849.0- 
400910 particularly interesting. Present data show that ab- 
sorption features in XDINs arc in all cases dependent on 
the spin phase. Until recently XDINs were believed to be 
steady (albeit some of them are pulsating) soft X-ray emit- 
ters. However, the brighter pulsating member of this class, 
RX J0720. 4-3125, has been discovered to undergo rather 
conspicuous spectral changes over a timescale of a few years 
(De Vries at al. 2004; Vink et al. 2004; Haberl et al. 2004). 
Weather this spectral evolution is accompanied by signifi- 
cant changes in the absorption line properties is not com- 
pletely clear as yet. 

Perhaps the most convincing evidence of long-term vari- 
ations of a cyclotron line comes from the SGR 1806-20. This 
source provided the first possible detection of an absorp- 
tion feature in a magnetar candidate (Ibrahim ct al. 2002; 
Ibrahim, Swank & Parke 2003). The feature has an equiva- 
lent width EW~ 0.5 keV, width cr ~ 0.2 keV and an energy 
of ~ 5 keV. The spectral feature appears to be transient 



4 INTERPRETATION IN THE MAGNETAR 
SCENARIO 

We discuss here the possibility that the onset of a twist 
in the magnctosphcrc might have been responsible for the 
two successive glitches, for the observed correlation between 
spectral hardening, luminosity and glitching activity, and, 
possibly, also for the transient behavior of the cyclotron fea/- 
ture. 

Thompson, Lyutikov & Kulkarni (2002) recently pro- 
posed a scenario in which the magnetars (AXPs and SGRs) 
differ from standard radio pulsar since their magnetic field 
is globally twisted inside the star, up to a strength of about 
10 times the external dipole and, at intervals, it can twist up 
the external fleld. The resulting stresses building up in the 
neutron star crust and the magnetic footprints movements 
can lead to crustal fractures, known to be responsible of 
the glitching activity of the pulsars (Rudcrman 1991b; Ru- 
dcrnian, Zhu & Chen, 1998, Ruderman 2004)^. These frac- 
tures could in principle be violent enough to affect substan- 
tially the neutron star interior and in particular to unpin 
the crustal neutron superfiuid vortex lines from the crustal 
lattice, causing the occurrence of a glitch (see Thompson & 
Duncan 1993). 

In this picture the vortex lines become unpinned se- 
quentially in small patches of the crust rather than in the 
whole crust as in the case proposed for radio pulsars (Ru- 
dcrman 1991b). Because of the large discontinuities of the 
crustal field, the reconnection of the vortex unpinned lines 
could be suppressed at such discontinuities, driving a pure 
horizontal field evolution to a point where reconnection is 
possible. This, in turn, can cause a second and potentially 
more energetic rupture of the crust (Thompson & Dun- 
can 1993). This scenario, therefore, can qualitatively explain 
the occurrence of the two successive glitches observed from 
IRXS J170849.0-400910 between 1999 and 2001 (Kaspi, 
Lackey & Chakrabarty 2001; DaU'Osso et al. 2003; Kaspi 
& Gavrifl 2003). We note that IRXS J170849.0-400910 was 
not totally recovered from the second glitch at the time of 
the BeppoSAX observation (Rea et al. 2003), while it was 
at time of the XMM-Newton observation reported here. 

More interestingly, the onset of a twist in the magne- 
tosphere can in principle explain both the observed spectral 



The strcngtli of tlic crustal luagnctic field can exceed the critical 
strength at which anisotropic magnetic field stresses break the 
crust. The evolving magnetic field is exijccted to fracture the crust 
at intervals, in a region of order of the coherence field length (for 
details see Flowers & Ruderman 1977; Romani 1990; Thompson 
& Duncan 1993). 
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Figure 6. Correlated photon index and X-ray flux changes with time. Dashed line represents the average source flux measured by 
previous ASCA (Israel et al. 1999) and ROSAT (Sugizaki et al. 1997) observations. The two arrows indicate the two glitch epochs. All 
fluxes are unabsorbed and calculated in the 0.5-10 keV energy band. The Chandra-HETG point (here we report the first order) comes 
from the analysis of an unpublished observation available in the Chandra archive (ObsID 2757). 



softening and a transient appearance of the cyclotron fea- 
ture. 

In fact, a key feature of twisted, force- free magneto- 
spheres is that they support current flows. The presence of 
charged particle (electrons and ions) produce both a large 
resonant scattering depth and an extra heating of the star 
surface, the latter induced by returning currents. Because 
the electrons distribution is spatially extended and the res- 
onant frequency depends on the local value of B, repeated 
scatterings could lead to the formation of a high-energy tail 
instead of a narrow line at the cyclotron frequency. In this 
model, both the scattering depth (which is strongly depen- 
dent on the magnetic latitude) and the luminosity released 
by particles hitting the surface increase with the twist an- 
gle. Therefore, since the spectral hardness increases with the 
depth, this implies a positive correlation between the source 
X-ray luminosity and the spectral hardening. This is indeed 
what the data reported in Fig|S| indicate. 

Moreover, the same model can provide a possible ex- 
planation for the transient appearance of a proton cyclotron 
feature during the epoch in which the twist was substantial. 
In fact, the charges present in the magnetosphere are also 
providing a large optical depth to resonant proton (ions) 
cyclotron scattering and the proton resonance occurs much 
closer to the star surface than the electron resonance. A tran- 
sient absorption feature at the proton cyclotron frequency 
may arise then when a) the twist angle is large and b) the 
X-ray luminosity at the resonant frequency (evaluated at 
the star surface, ujb,p) is large enough to exceed the gravity 
pull on the protons, so that positive charges are efficiently 
confined in a thin layer close to the star surface (see again 
Thompson, Lyutikov & Kulkarni 2002). 

Although not strictly compelling, the correlation be- 
tween line appearance, glitching activity and flux enhance- 
ment is therefore intriguing and suggests that the conditions 



for line formation may have been met at the epoch of the 
BeppoSAX pointing. 



5 CONCLUSION 

Our present analysis of the XMM-Newton observation per- 
formed on 2003 August of IRXS J170849.0-400910 revealed 
that this source is variable in X-ray over a period of about 
five years. This provides further evidence that AXPs, which 
were believed to be steady sources for a long time, exhibit 
both fiux and spectral changes. We pointed out the exis- 
tence of a clear correlation between the source flux level and 
spectral hardening. An interesting possibility is that both 
the increased X-ray emission and the hardening are related 
to the occurrence of glitches. 

The XMM-Newton pointing caught the source at a 
factor of two lower flux than that of the BeppoSAX ob- 
servation, performed soon after the last glitch, and with 
a significantly softer spectrum. By comparing the evolu- 
tion of IRXS J170849.0-400910 with that of the AXP 
IE 2259+586 , which followed a similar pattern, it is possible 
that IRXS J170849.0-400910 too had a transient period of 
bursting activity, missed in the RXTE monitoring. 

We interpret both the occurrence of the glitches and the 
spectral softening accompanying the flux decay in the frame- 
work of the magnetar scenario, by the onset of a twist in the 
external magnetic field (Thompson, Lyutikov & Kulkarni 
2002). Moreover, the occurrence of the twist could also ac- 
count for the possible variability of the cyclotron line, which 
was detected only while the source had not entirely recov- 
ered from the last, and more powerful glitch event. 
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